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Abstract: In heavy-ion collisions, quark-gluon plasma is likely to be produced with sizable 
initial pressure anisotropy, which may leave an imprint on electromagnetic observables. In 
order to model a strongly coupled anisotropic plasma, we use the AdS/CFT correspondence 
to calculate the current-current correlator of a weakly gauged U(l) subgroup of R symmetry 
in an M = 4 super- Yang-Mills plasma with a (temporarily) fixed anisotropy. The dual 
geometry, obtained previously by Janik and Witaszczyk, contains a naked singularity which 
however permits purely infalling boundary conditions and therefore the usual definition of 
a retarded correlator. We obtain numerical results for the cases of wave vector parallel and 
orthogonal to the direction of anisotropy, and we compare with previous isotropic results. 
In the (unphysical) limit of vanishing frequency (infinite time) we obtain a vanishing DC 
conductivity for any amount of anisotropy, but the anisotropic AC conductivities smoothly 
approach the isotropic case in the limit of high frequencies. We also discuss hard photon 
production from an anisotropic plasma and compare with existing hard-loop resummed 
calculations. 
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1. Introduction 

Information on the earliest stages of heavy-ion collisions and the production of quark- 
gluon plasma may be deduced from spectra of photons and dileptons, because real or 
virtual photons interact only weakly with the surrounding strongly coupled matter ||. |2|] . 
Thermal photon and dilepton production from an isotropic plasma has been studied both 
in (resummed) perturbation theory ||, H || and in strongly coupled super symmetric 
Yang- Mills plasma by means of gauge/gravity duality 0] (see also ||, |l0tl). 

The angular dependence of photon and dilepton production should carry information 
on anisotropies present in the plasma prior to isotropization [11, [l^, 13, 14, [15|, [16|, O], the 



details of which still remain to be understood. The success of fits of ideal-hydrodynamic 
simulations |l8|] suggested time scales of isotropization of < 0.7 fm/c, but later analysis 
using viscous hydrodynamics found that thermalization may be as late as ~ 2 fm/c, de- 
pending significantly on assumptions about the initial conditions [19]. The recent studies 
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of Ref. (2^, 21 1 suggest that large momentum-space anisotropics may be present for the 
most part of the time evolution both at weak and strong coupling. 

At parametrically weak coupling, a possible framework for carrying out systematic 
calculations in the nonequilibrium case of an initial momentum space anisotropy is given 
by the hard-loop effective theory |2^] where the backreaction of collective phenomena 
on the distribution of hard particles is ignored. The latter is either taken as station- 
ary and anisotropic |23|, |24], |25|, |26|, [27], EM EjJ, or free streaming with evolving degree of 



anisotropy J30|, 31, |32]|. Photon production from Compton-like and annihilation processes 



in a stationary anisotropic quark-gluon plasma has been studied in pi] ], confirming ex- 
pectations of a strong angular dependence. Combined with the time dependence of the 
momentum space anisotropy, this may even lead to yoctosecond pulses of high-energy 
photons f33]|. Weak-coupling results also exist for energetic dilepton emission from an 



anisotropic plasma [14], however without inclusion of next-to- leading order effects that 
become important when dilepton energies are not much higher than temperatures. 

In this paper we shall use the AdS / CFT correspondence to study electromagnetic prop- 
erties of a strongly coupled anisotropic plasma in the approximation of a fixed anisotropy. 
Like in the weak-coupling situation this should be a reasonable approximation for pro- 
cesses on time-scales smaller than the isotropization time, such as the production of hard 
photons or dilepton pairs. Following Ref. [34] we study the dual geometry of a stationary 
anisotropic N = 4 supersymmetric Yang-Mills plasma, which contains a naked anisotropic 
singularity that however still permits purely infalling boundary conditions. We use this 
setup to generalize the calculation of U(l) current-current spectral functions of Ref. 0] to 
the anisotropic case. Unlike Ref. ]^4| we refrain from expansions in the anisotropy param- 
eter, since even the smallest anisotropy changes the character of the fluctuation equations. 
Indeed, the spectral functions are strongly modified at small frequencies such that the DC 
conductivity vanishes. The absence of a horizon at any nonzero anisotropy leads to the 
absence of hydrodynamic behavior, which however just reflects the fact that a stationary 
anisotropy can only be a good approximation of the nonequilibrium dynamics at suffi- 
ciently short time scales. In fact, at large frequencies our results for the spectral function 
smoothly approach the isotropic limit. This suggests that our setup should permit us to 
study the effects of temporary anisotropics in a strongly coupled plasma on electromagnetic 
observables such as hard photon and dilepton spectra. 



2. The dual geometry 



There has been much progress recently in studying, numerically, the gravity dual of the 
nonequilibrium dynamics corresponding to analogs of heavy-ion collisions and quark-gluon 
plasma formation p! 



36 



39 



41, 42 



building upon and generalizing the 
AdS/CFT correspondence of maximally supersymmetric Yang-Mills theory at infinite 't 
Hooft coupling Ej|. Even though all these studies indicate very quick transition to hy- 



drodynamic behavior, estimated to take place at ~ 0.3 fm/c, the results of Ref. [42] show 
significant anisotropics in the energy-momentum tensor in the subsequent evolution. This 
is qualitatively in agreement with the analysis of Ref. [ 20 , EJl | studying kinetic theory equa- 
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tions in the presence of momentum space anisotropics at various amounts of specific shear 
viscosity with the result that both at strong and weak coupling one can expect sizeable 
pressure anisotropies. Naturally, at weaker coupling (larger viscosity), the anisotropy of 
the energy-momentum tensor reaches larger degrees and persists for longer times. 

Since the AdS/CFT study of electromagnetic spectral functions may be prohibitively 
difficult in a time-dependent gravity background involving the formation of (time-dependent) 
horizons, we shall content ourselves with a stationary anisotropic background obtained pre- 
viously by Janik and Witaszczyk [34]. 



2.1 Review of the Janik- Witaszczyk solution 

According to the AdS/CFT dictionary, the metric of asymptotically anti-de Sitter space 
in Fefferman-Graham coordinates 

dg2 _ ^f,u(x a ,z)dx^dx u + dz 2 
z 2 

where z is the holographic coordinate, encodes the energy-momentum tensor {T fJiU (x a )) of 
the dual super- Yang-Mills theory through its asymptotic behavior near the boundary at 
z = 0, 

7 ^(x CT , z) = + A$CO + 0(z 6 ) (2.2) 

with 



<^(* CT )> = StIIV)- (2-3) 



N 2 
2^ 2 

This fixes the boundary condition for the 5-dimensional Einstein equations with a negative 
cosmological constant 

Rmn(x, z) = -Ag M N{x, z) (2.4) 

and one can solve for the 5-dimensional metric. Requiring nonsingularity of the resulting ge- 
ometry would allow one to select physical solutions for the profile of the energy-momentum 



tensor and its time evolution p5l 44, 45 1 



However we are interested in a stationary anisotropic energy momentum tensor 

{TfMu) = diag(e, P L , P T ,P T ), (T%) = 0, (2.5) 

which we consider as an approximation to the full dynamics at sufficiently short time scales. 
With these boundary conditions the solution for the metric takes the form 



ds 2 = ( -a(z)dt 2 + b(z)dx\ + c{z)ds^ + dz 2 ). (2.6) 



with 



a[z 



) = (J _|_ ^2 z 4y/2-V36-2B 2 /4 _ ^2^1/2+^36-2574 



= (1 + A 2 Z 4 ) 1/2 ~ B/3+v/36 ~ 2B2 /12 (I - ^2^4^1/2+5/3-^36-252/12 ^ 7) 

c f z \ = (1 + ^2^1/2+5/6+^36-252/12^ _ ^2^1/2-5/6-^36-252/12^ 
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Table 1: Relation between different anisotropy parameters 
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The parameters A and B are related to the energy density and the pressures according 



to 
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2A 2 B 
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(2.8) 



The anisotropy of the system is parametrized by the dimensionless parameter B. For 
B = we obtain the usual anti-de Sitter/Schwarzschild solution in Feffer man- Graham 
coordinates. For positive (negative) values of B we obtain an oblate (prolate) anisotropy 
of the system. Pl and Pt are positive quantities in the range -V6 < B < V2, while in the 
larger interval — \/l8 < B < vl8 one of the pressure components can be made arbitrarily 
large and negative. The dimensionful parameter A, which in the isotropic limit equals 
7T 2 T 2 /2, will be set to 1 in the following discussion. 

In Table |l| we display the dependence of the pressure components normalized to the 
energy density as a function of B for a number of values used further below in numerical 
evaluations. Also given is the connection with the anisotropy parameter £ introduced in 
Ref. [46] through one-dimensional deformations of isotropic particle distributions /j so in 



momentum space, 



p 2 + £p| 



(2.c 



with jV(£) some normalization factor (e.g., M = y/1 + £ if the number density is to be kept 
fixed), and the corresponding energy- momentum tensor in kinetic theory 



T ».v _ N 



eft'. 



d 3 



p 



(2tt): 



■pV7(p)- 



(2.10) 



(This of course covers only the range where both Pl and Pt remain positive.) 
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Figure 1: Asymptotically spherical congruences of (holographically) radial light-like geodesies 
which get deformed into ellipsoids as they approach the singularity at u — 1 in units where A = 1. 
The blue (darker) surface corresponds to prolate anisotropy B = —yE, the red (lighter) surface to 
oblate anisotropy B = \/2, and the transparent mesh to the isotropic case B = 0. Here xt and 
xl correspond to the transverse and longitudinal spatial extents of the ellipsoids, which degenerate 
into an infinite disk or line for oblate or prolate anisotropy, respectively. (Note that xt is a radial 
variable.) 



Using the analytic expressions for the corresponding e, Pl and Pt of Ref. 
relation between £ and B is given by 



V36 - 2B 2 + 2B £ - 1 



+ 



V36 - 2B 2 - AB 2 ' (£ + l)£~i/2 atan^/2 _ i' 
For small anisotropics £ = %B + 0(B 2 ). 

2.2 The appearance of naked singularities 

For later convenience we write the metric in terms of u = z 2 : 

ds 2 = g tt {u)dt 2 + g L L(u)dx 2 L + g TT (u)dx T + g uu (u)du 2 , 



§1|, the 
(2.11) 



(2-12) 



gtt ( u ) = -1(1 + A 2 u 2 ) 1 l 2 - V ™ Z2 ~ m '\l- A 2 u 2 )y 2+VMZ2W / 4 
u 

g LL ( U ) = -(1 + y l2 lt 2Nl/2-B/3+V36-2BVl2 M _ ^2^2^1/2+5/3-^36-2^/12 
U 

gTT ^ = + A 2^2) 1/2+5/6+^36=25^/12 ^_ ^2^2^1/2-5/6-^36=2^/12 



U 



(2-13) 



.(«) 



1 

4u 2 ' 
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The geometry corresponding to this metric is pathological in the sense that a naked 
singularity appears whenever B does not vanish. For instance, the induced metric at 
constant t and u = l/A is degenerate, 

9LL9 2 TT oc (1 - A 2 )l 6 -»]/ 4 (2.14) 

Three-dimensional space degenerates at u = l/A into a two-dimensional sheet when B > 0, 
and into a one-dimensional line when B < 0. This is illustrated in Fig. [I] in terms of 
asymptotically (u — > 0) spherical congruences of holographically radial light-like geodesies. 
At finite u these are deformed into ellipsoids, which degenerate at u = l/A. 



However, in [34] it was noted that it is nevertheless possible to define purely ingoing 
and outgoing boundary conditions at the naked singularity, 1 which gives the possibility to 
define retarded correlation functions along the lines of Ref. [^7]] . In the following we shall 
do so for current-current correlators of a U(l) current selected from the -R-charge currents 
and assumed to be weakly gauged and thus coupled to electromagnetism. 

3. Current-current correlators and photon production 

Photons and leptons produced in an early stage after the heavy ion collision interact only 
very weakly with the medium and can leave the plasma without any rescattering. Therefore 
the production of photons and dileptons in a strongly coupled plasma is an interesting 
observable to study the far from equilibrium properties. In isotropic J\f = 4 SYM plasmas 
the emission spectra of photons and leptons at strong and weak coupling have been studied 
in M, In the following we shall closely follow this paper and generalize to the above 
anisotropic background. 

3.1 Basic facts 

The rate of photons produced per unit time and per unit volume is 2 

= 7^^^C<{K)\ k0=M , (3.1) 



where 



C< U (K) = J d±Xe- iK - x (J* M (0)J? M (X)) (3.2) 



is the Wightman function of electromagnetic currents. Any correlator of the conserved 
current can be expressed in terms of the spectral function, which is defined as 

XfW (K) = [ d*Xe- iK - x {[J* M (S>), J* M {X)]). (3.3) 



x But the expansion in small B performed in |}3] is not really allowed, because the character of the 
singularities in the equations of motion changes with the anisotropy parameter. We thank Karl Landsteiner 
for pointing this out to us. 

2 We use capital letter to denote 4 vectors and lower case letters for the absolute value of 3 vectors. The 
Minkowski metric is defined as r/^ = diag(— 1, 1, 1, 1). 
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Later we will determine this spectral function by making use of 

X»u(K) = -2ImC;f(K). (3.4) 

The Wightman function is related to the spectral function by 

C<(K) = -ix„u(K)(f>(K), (3.5) 

where <f> reduces to the Bose-Einstein distribution in thermal equilibrium. 

The dilepton production takes place via an intermediate virtual photon and the rate 
is given by 

where e; is the charge of the lepton. Here the Wightman function must be evaluated for 
time-like momenta. 

3.2 Introduction of photons and leptons 

After the brief recapitulation of photon and dilepton production rates in thermal media, we 
quickly review the possibilities to couple electromagnetism to the M = 4 SYM theory, which 
consists of SU{N C ) gauge bosons, four Weyl fermions ip p and six real scalars cp pq = —4> qp , 
all transforming in the adjoint representation of SU(N C ). There is also an anomaly free 
global SU(4) i?-symmetry present, under which the fermions transform in the 4 and the 
scalars in the 6. We can take a £7(1) subgroup of the ££7(4) i?-symmetry associated with 
a £7(1) gauge field coupled to the conserved current. By doing so we are able to model 
the electromagnetic interactions. In principle it is possible to take any linear combination 
of Cartan subalgebra generators to embed the £7(1) in the 5C/(4) i?-symmetry group (for 
details see J*]]). We will choose t 3 = (1/2,-1/2,0,0) to be the generator of the £7(1) such 
that two of the Weyl fermions have charge ±1/2 and two complex scalars have charge 1/2. 
The conserved current is 

p=3,4 

where a is the SU(N C ) group index and the covariant derivative D p involves the SU(N C ) 
gauge fields as well as the £7(1) electromagnetic vector potential A^. Because we are only 
interested in the leading order terms in the electromagnetic coupling e, it is sufficient to 
treat the electromagnetic interaction as being linear in A^. Then we can consistently ignore 
the electromagnetic vector potential in the covariant derivative acting on the scalar. In 
order to also add weakly coupled leptons I with charge e\ and mass m, the Lagrangian is 
extended to 

C = CsYM + C int --F^-l{p + m)l with £ int = eJ*A», (3.8) 

where J 3 is the i 3 component of the i?-current. 

In the gauge/gravity setup, only the SYM part will be realized dynamically, but we 
can calculate current-current correlators to leading order in the electromagnetic coupling 
and to all orders in the SYM coupling. 
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3.3 Tensor structure of anisotropic correlators 

For an anisotropic medium, the tensor structure of the current-current correlator is more 
complicated than in the isotropic finite temperature case. Because the electromagnetic 
current is conserved it must satisfy Ward identities, which implies that any Cau(K) ~ 
P^f{K) with P„ v = -q^ - K^K V /K 2 . 

We begin with the case when the wave vector is pointing in the direction of the 
anisotropy denoted by n = ne^. Then it is sufficient to introduce longitudinal and trans- 
verse projectors 

pT _ pT _ n pT _ r . _ P L - P P T C\ Q\ 

For the correlators we then find 

C^(K) = P^ T (K) + P^U L (K). (3.10) 

This is exactly the same structure one has in isotropic systems. 

If we choose the wave vector to point in a perpendicular direction with respect to the 
anisotropy k = ke^, we need one further tensorial structure, 

M)0 = "iO = 0, Pij = Oij p (3-11) 

MX) = MiO = 0' Pij = n 2 (3-12) 

Pfj,v = ~ P^v ~ Pfiw (3.13) 
As a consequence the correlator is then specified by three scalar functions 

C^{K) = Pl v U\K) + P* U I1 2 (K) + I%U L (K). (3.14) 

With a generic orientation of the wave vector, more structure functions would come 
into the play. We shall restrict our attention to the two extreme cases of wave vector 
parallel and orthogonal to the anisotropy direction. It is plausible that the generic case 
will interpolate smoothly between those two. 

3.4 Equations of motion and asymptotic solution 

To obtain the retarded correlator we have to solve the equations of motion of a gauge field 
in the geometry described by ( 2,12| ), which are given by 

d A (V=99 AC 9 BD F CD ) = 0. (3.15) 

At first we consider the case when the wave vector points in the direction of the anisotropy 
(k = &L e L)- Then the equation of motion for Er = ojAt is 

8 ( r^„uu TT\ tt 2 , n LL h 2 

E't + M ^± ° ] E' T - 9 ±S h E T = 0, (3.16) 



_ gg UUgTT 1 gUU 
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where primes denote derivatives with respect to u and the dependence of u is suppressed 
everywhere. For the longitudinal electric field El = k^At + ojAl we find 



' 1 ^g uu g tt g LL (g tt cu 2 + g LL k 2 L ) L 

g tt u 2 + g LL k 2 



L) 

tt, ,2 i „LL 1.2 

E L = 0. (3.17) 



gUU 



If the wave vector points in a direction perpendicular to the anisotropy, we obtain 
three different differential equations. Due to symmetry we are free to choose k = k±e± and 
then find 



2 V~99 UU 9 TT 2 9 UU 2 



(3.18) 



// d u (J^g uu g LL )^ g tt u 2 + g TT k 2 , N 

w/ u\y yy y j g _ y y 1 £ = Q /g ^ 

for the two modes transverse to the wave vector and 

F „ (g a ) 2 dn(^"V r V + (g rT ) 2 94^g""g a )fc 2 p 

i + V=99 uu 9 tt 9 TT (9 tt u 2 + ff TT A:f ) 1 

tt 2 , n TT k 2 

_ g u +g k, = 

gUU * V 7 

for the electric field along the direction of the wave vector. These equations are quite 
lengthy if the explicit form of the metric coefficients is inserted, therefore we will not do 
so here. 

For all of these equations, we can use a Frobenius ansatz near the boundary (u = 0) 
and find the characteristic exponents to be and 1. Close to the naked singularity (which 
appears exactly where the horizon is in the isotropic case) the differential equations have 
the following form, 

£* + c£os* + (r^*= (3 - 21) 



with a = (2 + y/36 — 2£? 2 )/4 < 2. For isotropic systems a = 2 and a Frobenius ansatz 
is still possible about u = 1. In that case we find the characteristic exponents ±ic*//\/8 
near the horizon and we can easily define ingoing boundary conditions. For nonvanishing 
anisotropy a < 2 and we can perform a coordinate transformation x = (1 — u)( 2 ~ a ^ in order 
to find appropriate boundary conditions at the naked singularity 3 . Then the equation of 
motion is given by 

d 2 ± p d ± 7 2 



+ -—<p + —<f> = (3.22) 



dx 2 x dx 



x 



where 



7 ~j2^af- (3 - 23) 



3 Thanks to Karl Landsteiner for clarifying this point. 
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The solution to this differential equation can be written as 

<j>(u) ~ (1 - n)( 2 - Q )( 1 -' 3 )/ 2 H; i 4 ) (2 7 (l - u)^ 2 ), (3.24) 

(2) 

where the Hankel function of the second kind HJ> represents ingoing boundary conditions, 
which we will use for our numerical studies later on. 

3.5 Spectral functions at strong coupling 

In order to find the on-shell boundary action, we start from the five-dimensional Maxwell 



action, which is given by 



S 5D = --^r I ^g~g AC g BD F AB F CD (3.25) 



4s? 



'B 

with g B = 16ir 2 R/N 2 and R the AdS radius. Choosing the gauge A u = we obtain the 
on-shell boundary term 



s * = -Jt I ^gg uu (g tt A' t {K } u)A t (-K,u) 
2g B Ju^o v 

+g LL A' L (K,u)A L (-K,u) + g TT A' T (K,u)A T (-K,u) ) (3.20) 



1 LL A' L (K, u)A L (-K, u) + g TT A' T 

Considering a wave vector pointing into the direction of the anisotropy first, we obtain the 
relation 

ug tt A' t (K, u) - kg LL A' L (K, u) = 0, (3.27) 

which follows directly from the equations of motion. Inserting this into ( 3.26| ) and rewriting 
the result in terms of electric fields El = ujAl + kAt and = ojAt, the boundary term 
of the action becomes 

TT \ 

+^-B' T (K, u)B T (-K, u) . (3.28) 

The transverse correlator is defined as 

5 2 S B _ u 2 5 2 S B 

Tt( ' ~ 5A T (K)5A T (-K) ~ 5E T (K)5E T (-K)' (3 ' 29) 



Applying the Lorentzian AdS/CFT prescription [47]] and inserting the explicit form of the 
metric coefficients we find 



2 E'(K,u) 
— lim — 

g B u^o Et(K,u) 



Ctt(K) = U T t(K) = — T lim J_ ( J _ [ ■ (3.30) 



After a similar computation for the longitudinal correlator we find 



U L (K) = -A Hm |^4- ( 3 - 31 ) 

g B u ^°E L (K,u) 
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Figure 2: Transverse contribution to spectral density for light-like momenta with values for the 
anisotropy parameter B — (black), B = 0.1 (red, dashed), B = —0.1 (red, dotted), B = 1 (blue, 
dashed), B = -1 (blue, dotted), B = V? ( green, dashed) and B = — V6 (green, dotted). The 
dimensionful parameter A, which equals 7r 2 T 2 /2 in the isotropic case B = 0, has been set to unity. 



When the wave vector is in the 1-direction we obtain three scalar functions which are 
given by 



2 , E'JK.u) 
U a (K) = — =- lim 



(3.32) 



with a = 1, 2, L. Notice that the equations of motion for El differ for wave vector perpen- 
dicular and parallel to the direction of anisotropy, cf. (3.17) and (3.19|). 



4. Numerical results 

4.1 Wave vector parallel to anisotropy direction 

First we discuss the form of the spectral function for a wave vector parallel to the anisotropy 
direction. When we consider light-like momenta only the transverse contribution to the 
spectral density is nonvanishing. The results for different anisotropy parameters are shown 
in Fig. ||. The black line corresponds to the isotropic case and coincides with the spec- 
tral density presented in after the correct normalization is chosen. For nonvanishing 
anisotropy we notice a qualitative difference for small frequencies, namely all spectral func- 
tions tend to zero faster than u, which is in stark contrast to the isotropic situation. How- 
ever for small values of B the curve quickly rises and then settles close to the isotropic one. 
For increasing anisotropy the spectral densities deviate more strongly from the isotropic 
result, in particular at small frequencies. 

Since small frequencies correspond to long time scales, the approximation of a sta- 
tionary anisotropy is bound to break down in this limit. On the other hand, at larger 
frequencies our analysis should be able to map the nonequilibrium situation in the form 
of a snapshot, so it is reassuring that there the effects of the anisotropy connect smoothly 
with the isotropic limit. Given that the appearance of the naked singularity is not a small 
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Figure 3: Double logarithmic plot to estimate the deviation from the linear behavior of the spectral 
function with respect tow. In the left panel we consider positive anisotropics and in the right panel 
negative ones. Color coding as in Fig. || 




modification as it changes the character of the differential equations, this is not completely 
obvious a priori. 

Since the low frequency limit seems to be unphysical, the question arises down to what 
value of oj we might trust this calculation. It makes sense to assume that the lower bound 
of the frequency will depend on the anisotropy of the system. For larger anisotropics our 
assumption of a time invariant background may break down earlier. A possible hint to 
estimate the range of validity of the calculation presented here can come from considering 
the deviation from x oc u; as is done in Fig. ||. We see that even for the most extreme 
anisotropy parameters this suggests our calculation to be valid down to w ~ 1 (in units 
where A = 1). 

In Fig. |] the form of the transverse part of the spectral function is shown for fixed 
and uj fixed, respectively. In the left panel the difference in the small frequency behavior 
between isotropic and anisotropic results is again obvious. 

As already mentioned the longitudinal part of spectral density vanishes for light-like 
momenta. This is still the case if we turn on the anisotropy. In Fig. [5] we see that the 
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spectral function is negative for space-like momenta, vanishes for light-like momenta and 
becomes positive for time-like momenta. 4 

4.2 Wave vector perpendicular to anisotropy direction 

Next we consider a wave vector pointing in the 1-direction. Then there is a mode which 
is both, perpendicular to the wave vector and the anisotropy direction. This mode can be 
compared to the transverse mode before and we see that the behavior for small frequencies 
and wave vectors is quite similar. The spectral function is larger for oblate anisotropy and 
smaller for prolate. However, while this is still true for larger frequencies and momenta 
in the present case (see Fig. ||] and the opposite was true for xt, where the behavior 
changed at some intermediate frequency or momentum. 

For k = k\e\ there is another mode which is transverse with respect to the wave vector. 
The results for this mode, which is pointing along the anisotropy direction are shown in 
Fig. |8] and H. Compared to the previous modes that were transverse with respect to the 
wave vector, the dependence on B changed. Here the spectral density is larger for prolate 
anisotropy and smaller in the oblate case. 

Finally we can consider the mode longitudinal to the wave vector. The behavior is 
shown in Fig. [[(]. The spectral density for light-like momenta vanishes, as it should. 

4.3 Conductivities 

The fact that for any nonzero anisotropy parameter B all spectral functions tend to zero 
stronger than linearly in the limit oj — > means that both the diffusion constant and 
the DC conductivity vanish according to Kubo's formulae. This absence of hydrodynamic 

4 It turned out to be hard to solve the differential equations in the longitudinal cases numerically. While 
the default setting of Mathematica's NDSolve produces error messages at some data points, but matches 
with most of the other methods whenever no error occurs, the so called "StiffnessSwitching" method pro- 
duces no errors at all but sometimes produces very different results from the other methods. We circum- 
vented the difficulties by using the default method and computing many data points, such that we can safely 
neglect all the data points producing an error message. The final plots do not change when the number of 
data points is further increased. 
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behavior is clearly related to the absence of a horizon, which for B ^ gets replaced by a 
naked singularity. 

In Fig. |ll] we display the results for the AC conductivities, juxtaposed for the cases 
of prolate (B < 0) and oblate (B > 0) anisotropies. In each case one can define longitudi- 
nal and transverse conductivities with respect to the direction of anisotropy. For prolate 
anisotropies transverse conductivities are found to be reduced compared to the isotropic 
case, whereas for oblate anisotropies this is true for longitudinal conductivities. However, 
in the limit of vanishing frequency, all conductivities go to zero. The frequency range in 
which this happens increases as the amount of anisotropy is increased, up to the point where 
one of the pressure components goes to zero. Curiously, when increasing the anisotropy 
parameter such that also negative pressures are produced, this trend is eventually reversed. 

4.4 Anisotropy of traced spectral function for light-like and time-like momenta 

The production of real photons is proportional to the trace of the spectral functions for 
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light-like momenta. When the wave vector is parallel to the anisotropy direction, we have 

X%(k = k L e L , K 2 = 0) = 2 XT (k L ), (4.1) 
whereas for transverse wave vector, 

X %(k = fciei, K 2 = 0) = X2 (fci) + Xl^x). (4.2) 



In Fig. 12 the results for the wave vector pointing parallel and perpendicular to the 



anisotropy direction are juxtaposed, and in Fig. 13 the ratio between the latter and the 
former is shown as a function of frequency and anisotropy parameter B. 

For frequencies oj < 1 (in units where A = 1) this ratio shows a rather dramatic de- 
pendence on the anisotropy parameter. However, as we have discussed above, we consider 
this regime to be unphysical since the assumption of stationarity of the anisotropic plasma 
ceases to make sense. For larger frequencies we indeed find a smoother dependence on the 
anisotropy parameter. For positive B (oblate anisotropy) the spectral function with wave 
vector pointing in the direction of anisotropy is reduced, for negative B the situation is 
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Figure 11: Anisotropic AC conductivities: prolate vs. oblate anisotropies for various values of B. 
Full lines correspond to longitudinal conductivity, dashed lines to transverse conductivity. Color 
coding of B as in Fig. ^ with the addition of orange lines for B = =f3, values for which negative 
pressures arise. 



reversed. As we shall discuss below, this is in line with the behavior of particle distribu- 
tions for corresponding momentum anisotropies, but the latter typically have exponential 
suppression at high momentum. 

In Fig. |l4] we display the behavior of x^n f° r time-like momenta, which is relevant 
for dilepton production, as a function of ki and k\ for the two cases of longitudinal and 
transverse momentum at a fixed value of K 2 = 1. 



4.5 Photon and dilepton emission 

In order to obtain the photon and dilepton production rates, we must insert the results 



for the spectral function into (3.1) and ( p.6[ ), respectively. For the final result we would 
also need to know the distribution function cf>, which in the out-of-equilibrium situation is 
no longer fixed to the Bose-Einstein function. Absent the fluctuation-dissipation theorem, 
we should therefore calculate the Wightman function within the gauge/gravity duality 
framework directly. Attempts to incorporate the full formalism of nonequilibrium physics 




Figure 12: X^n/u for light-like momenta with wave vector parallel to the anisotropy direction 
(left panel) and transverse to it (right panel). Color coding of B as in Fig. [ll] 




Figure 13: Ratio of x% for light- like momenta with transverse wave vector over this quantity 
with parallel wave vector (w.r.t. the direction of anisotropy) as a function of both frequency and 
anisotropy parameter B. 



in the AdS/CFT correspondence were undertaken in [pj| , [49| ], but we have not been able 
to apply these concepts to our case. 

We shall instead make an estimate of photon and dilepton emission rates by assuming 
that the distribution function (j) is given by the form (|2.9| ) with a parameter £ that at weak 
coupling gives the same energy momentum tensor as in the boundary field theory of our 



gravity duality (see eq. ( 2.11 )). Since the distribution function (|2.9|) depends exponentially 
on the angle, the overall normalization constant A/"(£) does not play an important role. 
Because the energy density at strong coupling, eq. ( |2.8|) , depends rather weakly on the 
anisotropy parameter B, we have fixed AA(£) such that the energy density at weak coupling 
remains constant as £ and B are varied. 

In Fig. 15 we display the resulting photon production rates in longitudinal and trans- 
verse directions for small and medium anisotropies. The left panel shows the situation for 



-17- 



5.0 



5.0 



4.5 ; 

4.0 ; 
3.5 ; 



4.5 - 
4.0 ; 
3.5 ; 



Figure 14: x M M for time-like momenta with invariant mass K 2 = 1 as a function of and k±, 
respectively. Color coding as in Fig. || 
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Figure 15: Photon production rate in forward (purple) and transverse (cyan) direction for B = 
±0.1 (dashed) and B = ±1 (dotted), where positive signs (oblate cases) correspond to the plot on 
the left and negative signs (prolate cases) are shown in the right panel. As a reference the isotropic 
rate is given in black. 



oblate anisotropics (B = 0.1 and 1), the right panel for prolate anisotropies (B = —0.1 and 
— 1). This agrees qualitatively with the corresponding result for oblate anisotropies in the 
weak coupling (hard-loop resummed) calculation of Ref. [jllj. 5 Similar results are obtained 
for production rates of dileptons. The angular dependence in both cases is dominated by 
the function <j) when it is chosen in accordance to the situation at weak coupling. 

5. Conclusion 

In this work we have generalized the AdS/CFT calculation of electromagnetic observables 
in strongly coupled Af = 4 super- Yang-Mills theory in thermal equilibrium of Ref. @ 

5 Our estimated result for real photon production is in fact more strongly suppressed in the forward 
direction than the result of Ref. Jllj]. It shares this behavior with the soft part of the full result of Ref. [O, 
which in the forward direction is dominated by hard contributions. The final result of Ref. [|ll| actually 
depends on the choice of a separation parameter of hard and soft scales which was fixed by a minimization 
procedure in the isotropic case. Fixing it instead anew for each value of £ would further reduce the hard-loop 
result in the forward direction. 
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to the case of a stationary anisotropic plasma using the singular geometry of Janik and 



Witaszczyk [34]. In contrast to the latter work, we have refrained from an expansion in 
the anisotropy parameter since any nonzero anisotropy changes qualitatively the character 
of the fluctuation equations. In the limit of small frequencies, we have found that for any 
amount of anisotropy the hydrodynamic behavior of spectral functions is lost, whereas at 
higher frequencies deviations from the isotropic case are smooth. This reflects the fact 
that a stationary anisotropic geometry can only be viewed as a reasonable approximation 
to the early stage of plasma formation and thermalization when the involved time scales 
are sufficiently short. 

We have worked out the independent spectral functions of the electromagnetic current- 
current correlation tensor (defined by a weakly gauged U(l) subgroup of the R symmetry) 
and their dependence on anisotropy parameter and wave vector. Absent a recipe for directly 
calculating Wightman functions, we have estimated photon production rates by adopting 
deformed thermal distributions that have previously been employed in weak coupling cal- 
culations. Doing so, our results turned out to be qualitatively similar to the weak coupling 
results of [11] in that an oblate anisotropy leads to a strong suppression in the forward 



direction. It would of course be very desirable to have a direct calculation of Wightman 
functions in our nonequilibrium geometry which so far has been achieved only in few cases 
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Recently a different gravity dual for anisotropic plasmas has been proposed in Refs. 



1 51, 52] which is based on geometries dual to Lifshitz-like fixed points constructed in Ref. 
1 53]. This approach involves additional bulk fields and thereby achieves completely regular 
geometries. It would be very interesting to compare spectral functions of current-current 
correlators obtained in this approach with our results. We would expect that the two 
approaches should agree in the regime of high frequencies but deviate in the infrared. 
However, in the latter regime neither of the two approaches will provide a consistent model 
of a strongly coupled anisotropic plasma, since the time dependence of the anisotropy can 
no longer be ignored. 
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